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The Planck and WMAP9 satellites, as well as the ATACAMA and South Pole telescopes, have recently 
presented results on the angular power spectrum of the comic microwave background. Data tentatively point to 
the existence of an extra radiation component in the early universe. Here, we show that this extra component can 
be mimicked by ordinary WIMP dark matter particles whose majority is cold, but with a small fraction being 
non-thermally produced in a relativistic state. We present a few example theories where this scenario is explicitly 
realized, and explore the relevant parameter space consistent with BBN, CMB and Structure Formation bounds. 
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I. INTRODUCTION 

The nature of dark matter (DM) is one of the major mys- 
teries in our current understanding of the universe. A com- 
prehensive strategy encompassing direct and indirect detec- 
tion plus collider searches, as well as indirect information 
from cosmological observations, will certainly be necessary 
to determine the fundamental nature of the DM. The re- 
sults from the DAMA/LIBRA Collaboration 1 1] and, more re- 
cently, from the CoGeNT |2, 3|, CRESST |4| and CDMS |5| 
experiments, might point to the first direct signals of light 
Weakly Interacting Massive Particles (WIMPs) scattering off 
of nuclei [6|. The gamma-ray excess from the Galactic Cen- 
ter region [7 1 and the tentative gamma-ray line observed with 
the Fermi Large Area Telescope fSl have also been interpreted 
as resulting from WIMP annihilation. The excess cosmic -ray 
positrons, recently confirmed by AMS-02 ||9], is an additional 
puzzling result that some have associated with WIMP annihi- 
lation [10-131 (see however |14|). 

Besides these putative direct and indirect detection signals, 
recent measurements of the angular power spectrum of the 
cosmic microwave background (CMB) by a variety of tele- 
scopes and satellites seem to indicate the existence of an extra 
component of radiation in the early universe. In particular, 
the WMAP Collaboration has recently presented their 9-year 
data and concluded that, after combining with data from Bary- 
onic Acoustic Oscillations (BAO), from the Atacama Cosmol- 
ogy Telescope (ACT), and from the Hubble Space Telescope 
(HST) (Ha = 73.8 ± 2.4 km s"^ Mpc"\ |15|), the number 
of effective massless neutrinos is A^e// = 3.84 ± 0.4 1.16.1 . 

The situation is, however, quite complex and subtle when 
the results of other experiments are considered. For exam- 
ple, the SPT collaboration has also recently reported new re- 
sults. Combining their data with BAO and HST, they con- 
cluded that Neff = 3.71 ± 0.35 |17|. Depending on the 
data set used, these conclusions can, however, qualitatively 
change. For instance, combining data from WMAP9 and SPT, 
Ref 1 18 1 finds Neff = 3.93±0.68, while from a combination 
of WMAP9, SPT, and the 3-year Supernova Legacy Survey 
(SNLS3) HU, Neff = 3.96 ± 0.69 ED. ACT also reported 
recently a much lower value than what quoted above, namely 
Neff = 2.79 ± 0.56 |20|, when utiHzing WMAP7 data, and 
Neff = 2.74 ± 0.47 when including WMAP9 data 1 18 1. This 
value can, however, be much larger if one includes BAO and 



HST data, reaching Neff = 3.43 ± 0.36 118|. It is clear that 
SPT and ACT results are somewhat in tension, and this may 
be related to the different values used for the lensing ampli- 
tude parameter Al |,21 J . 

Lastly, the eagerly awaited Planck results have just been re- 
ported. At face value, Planck data alone seem to disfavor an 
extra radiation component, pointing to Neff = 3.36 ± 0.34 
II22II . It is important to note, however, that the Planck col- 
laboration adopted fairly low values (at the 2.2a level) for 
Ho compared to previous studies |15|. Since Neff and Hq 
are positively coiTelated, increasing Hq would naturally yield 
higher values for Neff |22|. Interestingly, adding the Ho 
measurement and BAO data, the Planck Collaboration finds, 
in fact, Neff = 3. 52^0:45 E2- This latter value is - la 
above the Neff = 3.046 value favored by Big Bang Nucle- 
osynthesis (BBN) |23|. 

As the discussion above implies, the best fit value for Neff 
varies significantly depending on the data set used. Regard- 
less of that, the combined data reviewed above seem to be 
tentatively pointing towards the existence of an extra radia- 
tion component. Interestingly, it has been shown that this ad- 
ditional component can be interpreted as a thermally produced 
light DM species |24|, or as non-thermally produced relativis- 
tic DM particles, as shown in a model-independent way in 
Ref. 12511 . and in the context of supersymmetric frameworks 
inRef.l|26l. 

In this article, we present a few illustrative models, includ- 
ing a supersymmetric example, where the tentative extra ra- 
diation component is interpreted in terms of WIMP DM pri- 
marily composed of a standard cold component plus a small 
fraction (roughly 1% or less) produced non-thermally in a 
relativistic state. In Section III] we relate the energy density 
associated with non-thermally, relativistically produced DM 
with the number of effective neutrinos. In Section |lll] we ex- 
amine the cosmological bounds from BBN, structure forma- 
tion and the CMB on this scenario. In Section [TV] we present 
four models that satify the cosmological bounds and produce 
ANeff ~ 1. We summarize and conclude in SectionM 



II. LIGHT DARK MATTER MIMICKING AN EXTRA 
NEUTRINO IN THE EARLY UNIVERSE 

In order to relate the energy density associated with non- 
thermally, relativistically produced DM with the number of 



effective neutrinos, we start by calculating the ratio between 
their respective energy densities. Since the Cold Dark Mat- 
ter (CDM) and neutrino densities are given by pdm = 
Pc^Idmo.^^ and pi, = pc^^a^^N^/i, at Matter-Radiation 
Equality (MRE), the ratio between neutrino and DM energy 
density is 
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where il^ ~ 4.84 x 10"^, ftoM ~ 0.227, N„ is the number 
of neutrinos, and qeq ~ 3 x 10~^ is the scale factor at MRE. 

For Ni, = 1, we thus find that the energy density of one 
neutrino is ~^ 16% of the CDM density. As a result, if DM 
particles had a kinetic energy equivalent to •jdm — 1-16 at 
MRE, this fraction would produce the same effect as an extra 
neutrino species in the expansion of the universe at MRE, as 
already pointed out in Ref. |25 1. 

Assuming that some fraction / of the DM in the universe 
is produced via the decay X' — > DM + 7, we can relate the 
boost factor of the DM particle produced in the decay pro- 
cess (see AppendixlAJi with the number of additional effective 
neutrinos in the universe via 
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Eq. Q illustrates that this scenario posseses three free param- 
eters: 

(i) the X' lifetime for the decay process, 

(ii) the mass ratio Mjjm/Mx', and 

(iii) the fraction / of the DM density produced via the decay. 

With specific choices for the 3 parameters, one can easily re- 
produce AN^frr ~ 1 and therefore mimic the effect of an 
extra effective relativistic degree of freedom |25|. 

It is important to point out that in the present setup most of 
the DM would still be produced thermally, for example via the 
standard thermal freeze-out picture, with only a small frac- 
tion being produced non-thermally and with a large kinetic 
energy. We show in the next section that the fraction of DM 
produced relativistically must be small (/ <C 1) due to the 
fact this component directly affects the formation of structure 
and other cosmological bounds. For example, the electromag- 
netic energy released in the decay could have catastrophic im- 
pacts on the abundance of light elements as synthesized in the 
Big Bang Nucleosynthesis (BBN) framework, as well as on 
the spectrum of the Cosmic Microwave Background radiation 
(CMB). We analyze in detail structure formation, BBN and 
CMB bounds on / and on the X' lifetime in the following 
section. 



III. COSMOLOGICAL BOUNDS 
A. Structure Formation 

It is well known that the primordial universe had small 
inhomogeneities in the density distribution of matter, which 
evolved linearly at first, while eventually growing non-linearly 
to form the structures observed today. In the standard picture 
of structure formation, the dark matter played a crucial role in 
seeding the growth of structure at an acceptable rate. A crucial 
ingredient to this picture is that the DM particles must have 
had a negligible kinetic energy at the time of structure forma- 
tion. If not, the formation of structures at large scales would 
be slowed down because DM particles with large kinetic ener- 
gies would not cluster at sufficiently small scales due to free- 
streaming. This slowing down the growth of structures would 
wash out matter density fluctuations. It is therefore critical to 
investigate the suppression on the growth of structure caused 
by the fraction of DM which is non-thermally produced in 
a relativistic state via the aforementioned decay process, and 
to derive quantitative bounds on the fraction of DM particles 
produced in this manner 

It has been shown that, at scales above the free-streaming 
length, relativistic DM particles effectively behave like cold 
DM. At small scales, the matter fluctuations of cold DM parti- 
cles is governed by a linear equation which, during the matter- 
dominated stage of the Universe, follows the power-law be- 
havior l,27J 
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In the equation above, aoo(/) is the growth rate of the cold 
DM field, which during the matter-dominated epoch is given, 
to first order, by |271 
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where / is, as above, the fraction of the DM density produced 
non-thermally in a relativisitic state, and where / ^ 1. It 
is then possible to determine the suppression on small-scales 
caused by the non-thermal production of relativistic DM par- 
ticles by comparing the matter fluctuation given in Eqs.([3]l-(|4]| 
with the pure cold DM case at MRE. To first order, we find 
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As mentioned above, it is important to notice that Eq. (|5]l is 
applicable only in the matter-dominated regime, and a more 
accurate calculation considering the effects of, e.g., a cosmo- 
logical constant, is out of the scope of this work. 

After combining measurements of the amplitude of matter 
fluctuations in the Universe on scales of 8h^^ Mpc from the 
WMAP9 results |16| with Lyman-alpha forest data ||28l, the 
resulting limit is g > 0.95. This limit implies, from Eq. (J5]l, 
that / < 0.01. 

In summary, structure formation constraints require that 
only a small fraction (less than 1%) of the DM in the Universe 



might have been produced in a relativistic state from a non- 
thermal process such as the decay we consider in the present 
framework. 



B. Big Bang Nucleosynthesis 

As mentioned above, the lifetime and the energy released 
by the mother particle X' are also constrained by BBN 
bounds. This is because the energy released at a given time 
in the history of the Universe might induce electromagnetic 
showers that create and/or destroy light elements synthesized 
in the early universe [29 1. In this subsection we analyze the 
possible effects of this scenario on BBN. 

Given the overall quantitative success of BBN in the "stan- 
dard cosmological model", the first requirement for new 
physics is that it not drastically alter any of the light elemen- 
tal abundances such as those of D, "'He, ^Li. Energy releases 
which occur long after BBN may, in principle, spoil the suc- 
cessful BBN predictions for light elemental abundances. The 
energy of photons produced in late decays of X' would have 
been rapidly redistributed through scattering off background 
photons i'^^BG -^ e+e") as well as through inverse Compton 
scattering (e 7 — ?► e 7) ll29l[30l . As a result, the constraints we 
obtain from BBN are, to an excellent approximation, indepen- 
dent of the initial energy distribution of the injected photons 
and are only sensitive to the total energy released in the decay 
process ll29ll . 

In order to derive BBN limits on the fraction of relativistic, 
non-thermally produced DM via the decay X' -^ DM+'j, we 
therefore need to calculate the total electromagnetic energy 
released. Let Y = n/n'P^'^^ , where n is the number density 
of particles of a particular species and n^^^^ is the number 
density of CMB photons, given by 
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The total electromagnetic energy released from the X' decay 
is thus Eem = E^Ydm- If for each X' particle we have 
the production of a DM particle plus a photon, then Yx' = 
Ky = Ye)m,t = Ydm,o, where Yda/.t determines the number 
density of particles at a time equal to the lifetime of X', and 
Ydm,o is the number density of DM particles today. 

We thus find that the normalized number density of DM 
particles is given by 



Ydm — 
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This expression evaluates to 
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The factor / showed up in Eq. (jSll because we assume here 
that only a fraction of the DM of the universe is produced 
in the decay process, while the majority of it is produced in 
a non-relativistically by some other mechanism that does not 
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FIG. 1. The parameter space defined by the mother particle X' life- 
time (tx') and tlie fraction of relativistically produced DM (/) times 
the mother-to-daughter mass ratio Mx' /Mdm , and constraints from 
BBN and CMB. The shaded regions show BBN bounds on the non- 
thermal production of DM via the decay X' — > DAI + 7, cor- 
responding to an excess relativistic degrees of freedom AN^ff = 
0.2, 0.4, 0.6, 0.8, 1. The green curve represents the CMB bound (re- 
gions to the right of the curve are ruled out). We assume Mx' 3> 
Mdm- 



contribute any significant energy release at BBN (for example, 
via a standard thermal freeze-out process). 

Since the photon energy produced in the decay is. 
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the total electromagnetic energy released is given by 
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In the limit Adx' ^ Mjjm we can directly relate the to- 
tal energy release given in Eq. ( 10 1 with the quantity / x 
{Mx' /Mdm) as well as with AN^fj, as given in Eq. (J2|. 

With this in mind, we can use the limits on the light el- 
emental abundances derived in Ref. f29\. Since their re- 
sults were obtained as a function of the total energy released, 
we straightforwardly convert them in terms of the quantity 
/ X {Mx' /Mdm)- We show our results in Fig.fTl where the 
shaded regions show portions of parameter space ruled out by 
BBN, as calculated in Ref. |29| employing the "baryometer" 
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The straight lines in- 



dicate combinations of the quantity / x Mx' /Mdm and tx' 
producing the AN^ff as in the labels. 

Fig.fllshows that the BBN bounds are weak for early decays 
because at early times the universe is hot and the initial pho- 
ton spectrum is rapidly thermalized, leaving just a few high- 
energy photons able to modify the light elemental abundances. 
However, for lifetimes longer than 10"* s, BBN excludes most 
of the relevant parameter space. 



C. Cosmic Microwave Background 

Similar to the BBN constraints, bounds from measurements 
of the CMB spectrum are largely independent of the injected 
photon energy spectrum. These bounds also depend primarily 
only on the total energy release. The key effect of the ad- 
ditional energy injection in the form of photons is related to 
spectral distortions to the CMB black-body spectrum 1 3 1 1 . 

For times t < lO'^ s the processes of bremsstrahlung, i.e. 
eX —i' eX-f (where X is an ion), Compton scattering and 
double Compton scattering ej — > 677 quickly thermalize the 
injected photon energy 1 311 . For t > 10'^ s, however, the 
bremsstrahlung and double Compton processes become inef- 
ficient, and the photon spectrum relaxes to a Bose-Einstein 
distribution with a chemical potential (/i) different from zero. 
Limits on fi can then be used to constrain this additional in- 
jected energy component. 

The change in the CMB spectrum caused by the energy 
injection is constrained by precise measurements from cur- 
rent satellites and telescopes. Following the procedure used 
in Ref. |31|, we calculated the limit (shown as a solid green 
curve in Fig. [T]i imposed by CMB data on the non-thermal 
production of DM at a given late time. In particular, measure- 
ments of the CMB spectrum constrain any perturbation which 
might cause it to depart from an almost perfect blackbody 
spectrum. These limits are expressed in terms of deviations 
of the chemical potential fi from 0. Current bounds imply that 
/i < 9 X 10^'' 1 32 33 1 . This upper limit can then be translated 
into a bound on the amount of energy released during the de- 
cay, providing a constrain on the quantity / x {Mx' /Mjjm) 
for a given lifetime. 

Fig. [T] clearly shows that the CMB bounds become more 
restrictive than those from BBN for lifetimes tx' longer than 
^ 10^ s, ruling out almost all the relevant parameter space 
to explain AN^ff with relativistic non-thermally produced 
WIMPs. 



IV. FOUR ILLUSTRATIVE DARK MATTER MODELS 

To summarize, we found thus far that constraints from 
structure formation, BBN and CMB require (i) / < 1% and 
(ii) a lifetime for the mother particle X' shorter than 10'' s. 
For concreteness, we present here four example dark matter 
models, including a supersymmetric case, where DM particles 
could be partly produced non-thermally in association with a 
photon. We also show that these models have the potential 
to reproduce AN^ff ~ 1 while obeying the cosmological 
bounds described in the previous section. Unless otherwise 
noted, we shall use the symbol AIx' for the mass of the mother 
particle and Mjjm for the mass of the DM particle, while we 
use different symbols to indicate the associated fields. 

Before discussing specific models, we note that, in general, 
we can recast Eq. (|2| as 
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This relation, in view of the constraints obtained above, 
namely that r > 10^ s and that / < 0.01, implies 
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> 4 X 10^ ANeff. 



(12) 
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We thus find that the ratio between the mother and the DM 
particle, for appreciable values of AN^jf, must be at least of 
the order of 10^. We then conclude that to produce a large 
enough number of effective relativistic degrees of freedom, 
the mother particle and the DM particle cannot be nearly de- 
generate. In particular, this rules out scenarios where the two 
relevant particles are expected to be close in mass, such as, 
for example, in any universal extra dimensional setup (see e.g. 
Ref. Ol). 



A. Spin-0 particle 

In this first case, the mother particle X' is a heavy scalar {S, 
with mass Mx') which decays into a spin-1, stable DM par- 
ticle (i?p, with mass Mom) plus a photon. The relevant La- 
grangian interaction term is given by the effective dimension- 
5 operator 



Ceff = \B^{d,S)F^^\ 



(13) 



where F^^^ is the U{1) electromagnetic field tensor 

The resulting lifetime of the mother particle in this case is 
found tobeQ 
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with A the cutoff scale of the effective theory, as in Eq. ( 13 1. 
The upper limit on the lifetime we derived above, r < 
10^ s, can then be recast as a lower limit on the mother parti- 
cle mass, for a given cutoff scale A. This lower limit on the 
mother particle mass can then also be, in turn, converted to a 
lower limit on the DM particle mass Mom, once values for 
the parameters / and AN^jf are fixed. Using the expression 
for the lifetime, along with Eq. (|2]), we obtain the requirement 
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In Fig. l2]we plot the lower limit for the DM mass as a func- 
tion of the cutoff scale A, for / = 0.01 and r = 10^ s. To 
obtain DM masses in the 10 GeV mass range in this effec- 
tive model would require A ^ 10^^ GeV. Had we included 



a coupling constant g in Eq. ( 13 1, we would have concluded 
that for Mom ^ 10 GeV, if A is of the order of the GUT 
scale (A ^ 10^^ GeV) the required coupling would have been 
g ^ 10-3. 



' As explained above, we are well in the limit where the mother parti- 
cle is much heavier than the DM particle, see Eq. jl2^; we thus neglect 
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FIG. 2. Lower limits on the DM mass Mdm for a model in which a 
heavy scalar particle decays into a spin-1 DM particle plus a photon, 
for different AN^ff. The relevant effective operator for this model 
is given in Eq. (fTsl. We assume in this figure / = 0.01 and a lifetime 

T = 10* S. 



FIG. 3. Lower limits on the DM mass for a model in which a heavy 
spin-1 boson decays into a spin-1 DM particle plus a photon, for 
different ANeff. The relevant operator for this effective model is 



given in Eq. 1 16 1. We employed / = 0.01 and a lifetime r — 10* s. 



Note that the x-axis is multiplied by a factor 10 



B. Spin-1 particle 

In this case, a heavy spin-1 boson (W^, again with mass 
Mx') decays into a spin-1 DM particle (B^, of mass Mdm) 
plus a photon. A Lagrangian term which mediates this decay 
through a renormalizable operator is 



Ceff = gW^B^F^"^ 



(16) 



where, again, F'^" is the Standard Model electromagnetism 
field tensor. The lifetime of the mother particle in this case is 
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We can again use the upper limit on the lifetime t ^ 10^ s 
and Eq. (J2]i to set a lower limit on the DM mass, 



M, 



DM 



< 



2.93 • 10-36 GeV 
? 



Vl04 s) 



1/2 



.^^eff, 



(18) 

We can thus calculate the values for g and for the DM particle 
mass producing a given number of effective relativistic de- 
grees of freedom AN^ff, once the fraction of DM produced 
non-thermally (/ = 0.01) and the lifetime (t = 10"* s) are 
set to values compatible with BBN and Structure Formation 
bounds. We show our results in Fig. [3] 

The figure implies that in this particular scenario a DM par- 
ticle with a mass in the ^ 10 GeV range might account for a 
ANeff ^ 1 extra radiation component for g ~ lO"'^'^. 



C. Spin- 1/2 particle 

In this case a heavy fermion (V', with mass Mx') decays 
into a spin- 1/2 DM particle (x, mass Mum) plus a photon 
according to the effective dimension-5 operator 
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Ceff = ^{i> a^. x) F^'" + h-c. 
The lifetime of the mother particle in this scenario is 



4.14 X 10" 



-24 



A y /GeV\ 
M^J \M^) 



(19) 



(20) 



We can again use the upper limit on the lifetime to set a 
lower limit on the DM particle mass. 
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This case closely mirrors the spin-0 example of Sec. IV A the 
only difference being a numerical factor of 32 in the lifetime. 
This translates to moving the DM mass upper limit down by 
a factor of 324/'^ « 3. Therefore we conclude that, similarly 
to the spin-0 case, the non-thermal production of a ~^ 10 GeV 
DM particle can mimic one extra relativistic neutrino in the 
early universe for A ~ lO'^'^ GeV, as shown in Fig.0 



D. A Supersymmetric Example: tlie Bino/Gravitino System 

Neutralinos are the mass eigenstates resulting from a mix- 
ture of neutral B-ino, W-ino, and Higgs-inos. In many su- 
persymmetric models, including certain supergravity scenar- 
ios with universal soft supersymmetry breaking masses at the 



10^ 



10^ 



> 

Q 



10- 



10"' 





y": 


ANeFf = 0.2 








ANeFf = 0.4 

— - ANeFf = 0.6 
ANeFf = 0.8 


^a^ 


ANeFf = 1 ^ 


"/r^^i^ 


^^'<''^^ ^ 


^ ..^ ' 


'•.X" 


^.'<'<^ J 


^ ..^ ''"^^ 




^0'^'?5^ 


- 


^ ,^.'^<y^ 




^-<^^ 


- 


^^-^-'ij?^ 


_ 


^ y,''^'^ 


- 


^'<''^^ 


f = 0.01 


-^ '^y^ 




C'^ 


T = io''s : 


i^ 




1 1 1 1 1 1 ' ' ' 1 ' ' ' ' 1 


1 1 1 1 1 1 III 



10- 



10^^ 



10^' 



10^1 10^' 

A (GeV) 



10^ 



10^* 



10- 



> 
o 



10-' 



10-' 









■ 


— 


- ANeFF = 


0.2 


: 


— 


-ANeFF = 


0.4 


■ 


... 


- ANefF = 


0.6 


y' 


— 


- ANefF = 


0.8 


y^ 




-ANeFF = 


1 




F = 


= 0.01 










y' 


-' y' y^ 


y 


^ 


y 


' ^>^ 












1 


1 III 





10' 



10^ 



10^ 



10" 



10= 



10"^ 



T (sec) 



FIG. 4. Lower limits on the DM mass for a model in which a heavy 
fermion decays into a spin- 1/2 DM particle plus a photon, for differ- 
ent AN^ff . The relevant operator for this effective model is given in 



Eq. 1 19 1. We employed / = 0.01 and a lifetime r = 10* 



FIG. 5. Lower limit on the gravitino mass, for a supersymmetric 
model where a fraction / of gravitinos are non-thermally produced 
by the decay of pure Binos. The relevant lifetime is derived from 
Eq. |22}. As in the previous figures, we have used / — 0.01 and . 



grand unification scale, the lightest neutralino is an almost 
pure Bino. The lightest neutralino is also the next-to-lightest 
supersymemtric particle with the lightest supersymmetric par- 
ticle being the gravitino. Binos decay into a gravitino-photon 
final state via the interaction Lagrangian term lf35l 



C 



-G^iY^Yh^Bt^F^p, 



(22) 



where M^, ~ 2.4 x 10^* GeV is the reduced Planck mass. In 
models with low-scale supersymmetry breaking Mg ^ A/g,, 
which is exactly the limit needed to realize the scenario we are 
interested in. In this case, from Eq. ( [22] i we find a neutralino 
lifetime of 



t{B -> -fG) ~ 750 s 



IkeV 



IGeVy 



Mi 



(23) 



Similarly to previous sections, we can convert the bound on 
the lifetime into a lower limit on the DM mass. 



Mg^ 



I T \l/2 / 

(4MeV)(^) f 



/ 



VAiV, 



e// 



5/3 



(24) 



Setting, as above, / = 0.01 and r ~ 10"' s, we show in Fig.p] 
that a gravitino with mass M^ ~ Mdm in the 2 to 20 keV 
range mimics the effect of an extra neutrino while still obey- 
ing cosmological bounds. We note that such light gravitinos 
are marginally consistent with bounds from small scale struc- 
ture (for example Lyman-alpha forest data, see e.g. Ref. |36l). 
Also, assuming a gluino mass in the TeV range, gravitinos in 
the mass range needed here would require a reheating temper- 
ature Tfj close to the electro-weak scale to avoid over-closing 
the universe, since the thermal production of gravitinos in the 
early universe ifJTll implies 



Tn 



100 GeV 



0.2 



lkcV\ 
MG ) 



(25) 



Such a low reheating temperature would rule out certain sce- 
narios for the production of the baryon asymmetry in the uni- 
verse, such as leptogenesis, but is in general not phenomeno- 
logically implausible. 

Finally, we note that the constraints on the lifetime and 



Eq. (23 I imply a lower limit on the bino mass, which must 



be larger than about 1 GeV. This is, of course, perfectly com- 
patible with the supersymmetric models relevant here. 



V. CONCLUSIONS 

Recent data show tentative evidence for excess effective ra- 
diation degrees of freedom in the early universe. In this study, 
we analyzed scenarios where such extra radiation stems from 
a sub-dominant population of WIMP DM particles which is 
produced relativistically from the decay of a heavier particle. 
Such decay could affect structure formation, the cosmic mi- 
crowave background, and the production of light elements in 
the early universe. This framework is therefore constrained in 
a model-independent way by various bounds from cosmolog- 
ical observations. 

We showed that structure formation limits the fraction of 
the DM energy density that can be produced in the aforemen- 
tioned relativistic state to less than a percent, while CMB and 
BBN force the decay when the DM is produced to occur with 
a lifetime shorter than about 10"* s. Given these two key con- 
straints, we showed that a large enough contribution to the 
effective relativistic degrees of freedom ANeff is achieved 
only for a relatively large mass hierarchy (on the order of 10^) 
between the mother {X') and daughter (DM) particle masses 
in the decay (for example, X' — >■ DAI + 7). 

We presented a set of illustrative particle DM models that 
implement this scenario, including effective theories as well 
as a supersymmetric example. For each example, we derived 



the relevant combination of parameters leading to a significant 
ANeff, while being consistent with all cosmological bounds 
from structure formation, BBN and CMB. 

As a final remark, we note that in the context of the present 
framework it is in principle easy to accommodate both: 

(i) Mom — 10 GeV, i.e. in the mass region of partic- 
ular interest today from both direct detection ll38l and 
indirect detection |7|, and which has been extensively 
investigated in many particle physics models I39-4T1 . 
or 

(ii) assume that the DM particle mass is in the classic 
WIMP ^ 100 GeV range, which is also particularly 
interesting from the standpoint of possible indirect de- 
tection signals Il7l l42ti44l . and which has also been ex- 
haustively explored in particle physics models as well, 
see e.g. Ref |45-48|. 

For any DM particle mass, the key requirements that the 
present setup imposes are to have: 

1. a large enough mother particle mass Mx' compared to 
the dark matter mass, to produce a sizable AN^ff ac- 
cording to Eq. (|2]), 

2. a short enough lifetime for decay into the DM particle 
(r < lO'^ s) to evade BBN and CMB constraints, and 

3. a small enough fraction of relativistically produced par- 
ticles (/ ^ 0.01) to evade structure formation bounds. 



Substituting Eq. ( |A3[ ) into Eq. ( |Al[ l and taking Ex' 

we get 



M 



X' 



IDM = 



M 



X' 



M 



DM 



2M'^ 



(A4) 



2Mdm 
Taking into account the expansion of the universe and assum- 
ing that all decays happen at the decay lifetime t — t, we 
find 



IDM 



M 



X' 



2M 



'J- DM 



Mom 
2M'x 



- 1 



(A5) 



At the time of decay, the universe was dominated by radiation 
with a = (f/io)^/^, thus 



IDM 



1 



1/2 



M 



X' 



Mdm 
2M'x 



1 



(A6) 



2Mdm 
Therefore, at matter-radiation equality {t = tEq) we obtain 



7DAf -^ 1+7.8x10" 



VlOSsJ 



1/2 



/ M 



X' 



M. 



DM 



\2M, 



DM 



2M'x 



- 1 



(A7) 

Eq. ( A7 I gives us the boost factor of the DM particle as a func- 
tion of time and of the mother-to-daughter mass ratio, for the 
case where all decays happen when the universe was radiation 
dominated. 

Since the additional radiation density reads 



Pextra = ./ x Pdm[idm -1), 



(A8) 
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Appendix A: Boost Factor 

The boost factor of a DM particle produced in the decay 
X' — )■ DM + 7 can be found using energy and momentum 
conservation. 



Ex' — E^ + MuAI "fDM, 

Px' =E^ + Pdm- 
Assuming that X' decays at rest, we find 

1 



P, 



DM 



- '^^' - 2Mx' 



(M 



MIm)- 



(Al) 
(A2) 



(A3) 



where / is the fraction of DM non-thermally produced in 
the decays, we conclude that the number of effective rela- 
tivistic degrees of freedom at matter-radiation equality asso- 
ciated with the non-thermal production of DM is AN^fj — 
Pextra/piiy, whcrc pi^y is the number density of one neutrino 
species at the same epoch. Using Eq. (flj, we find 



AN.ff = 



filDM - 1) 

0.16 



(A9) 



Thus, substituting Eq. ( A6 1 into Eq. ( A9 1, we obtain 



ANeff ~ 4.87 X 10^' 
Mx' 



2M 



DM 



_i_ Mdm 
~2M\ 



1/2 



-1 



X /. (AlO) 



This expression indicates that if some fraction of the DM of 
the universe is produced non-thermally via the decay X' -> 
DM + 7, for reasonable values of the lifetime and mass ratio 
(Mx' /Mom), this production mechanism produces an effec- 
tive excess of relativistic degrees of freedom. 
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